Quantifying intermittent transport in cell cytoplasm 
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Active cellular transport is a fundamental mechanism for protein and vesicle delivery, cell cycle 
and molecular degradation. Viruses can hijack the transport system and use it to reach the nucleus. 
Most transport processes consist of intermittent dynamics, where the motion of a particle, such 
as a virus, alternates between pure Brownian and directed movement along microtubules. In this 
communication, we estimate the mean time for particle to attach to a microtubule network. This 
computation leads to a coarse grained equation of the intermittent motion in radial and cylindrical 
geometries. Finally, by using the degradation activity inside the cytoplasm, we obtain refined 
asymptotic estimations for the probability and the mean time a virus reaches a small nuclear pore. 
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Introduction. Cell transport, which may involve vesicles 
or proteins is essential for cellular function and home- 
ostasis. In general free diffusion in the cell cytoplasm 
is not efficient and many particles such as large viruses 
cannot pass the crowded cytoplasm [l[ without hijack- 
ing the complex cellular transport machinery and use 
molecular motors, such as dyneins, to travel along mi- 
crotubules (MTs) toward the nucleus. Both vesicular and 
viral motions alternate intermittently between periods of 
free diffusion and directed motion along MTs |2|. Such 
viral trajectories have been recently monitored by using 
new imaging techniques in vivo [HHj]. 

The switch nature of the motion, imposes a complex 
behavior of the particle trajectories which depends on the 
number and distribution of MTs, the rate of binding and 
unbinding and the diffusion constant of the free particle. 
Some physical properties, such as the mean velocity of 
trajectories has been obtained for the motion in domain 
made of parallel strips, in which a random particle has 
a deterministic motion on the stripes and pure diffusion 
outside [51- In case of a population of motors, at equi- 
librium between free diffusion and bound on MTs, the 
motor distribution has been studied in cylindrical and 
radial geometries in 0, 0] ; the authors estimate the for- 
ward binding rate using Brownian simulations in Q and 
experimentally in Q. 

We consider here a particle x(i), which can be de- 
scribed using the stochastic rule: 

/2Ddw for x (t) free 
<lx = < (1) 
V for x (t) bound 

where w is a standard Brownian motion, D the diffusion 
constant and V the velocity of the directed motion along 
MTs. 

In this communication, we compute the mean first pas- 
sage time of a single particle to a population of MTs. We 
thus provide an analytical expression of the forward bind- 
ing rate of a motor to MTs in both radial and cylindrical 
geometries. Using the analytical expression of the for- 



ward binding rate, we propose a coarse-grained descrip- 
tion of a switch dynamical motion of a particle, which 
can either be a virus, a vesicle or a molecular motor. 
This description, which is the main result of our paper, 
is a fundamental step to estimate the probability and the 
mean time to arrive at a small target. Moreover, using 
this description, we obtain the steady state distribution 
of virus in MTs network without resorting to the assump- 
tion of a two-state model @, [1] ■ 

We thus compute an effective steady state drift b(x) 
such that the particle motion |T]) can be coarse-grained 
by the stochastic equation: 



dx = h(x)dt 



(2) 



Using results derived in [T(j , equation ((2|) and the degra- 
dation activity in the cytoplasm due to protease or lyso- 
some, we obtain asymptotic estimates of the probability 
and the mean time for a virus to reach a nuclear pore. 
The problem of finding a small target is ubiquitous in 
cellular biology and recent theoretical studies [ll|, [l2j 
suggest that the geometrical organization of the medium 
play a fundamental role in this search process. 
Mathematical Modeling. We represent the cell cytoplasm 
as a bounded domain f2, whose boundary <9f2 consists of 
the external membrane dQ ex t and the nuclear envelope, 
both of which form a reflecting boundary dN r for the 
trajectories of ([2|), except for small nuclear pores dN a , 
where they are absorbed. The ratio of boundary surface 



areas satisfies e 



\ dN a 



<C I. We model the virus degra- 



dation activity in the cell cytoplasm as a steady state 
killing rate k(x) for the trajectories of so the sur- 
vival probability density function (SPDF) is the solution 
of the Fokker-Planck equation [l3j 



^=DAp-V-bp 
at 

P (x, 0) = pi (x) 
with the boundary conditions: 
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p(x, t) = on dN a and J(x, i).n x = on dN r U 9f2 ex t(3) 
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FIG. 1: The fundamental step is represented with a dotted 
line; a virus starts at a position xo, diffuses freely, binds to a 
MT over a distance d m and is then released at a final position 
Xf. The solid line represents a trajectory generated by the 
steady state equation ([2}. In the parenthesis, we point out 
the mean times for each portion of trajectories. 



where n x denotes the normal derivative at a boundary 
point x. The flux density vector J(x, t) is defined as 



J(x, t) = -L>Vp(x, t) + b(x)p(x, t). 



(4) 



The probability P/v and the mean time tn that a tra- 
jectory of ([SJ) reaches dN a are given by the small hole 
theory for two-dimensional domains Q and drifts b (x) — 
-V$(x), as 
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(5) 
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Hereafter we derive explicitly the steady state drift b(x) 
as a function of some geometrical and dynamical param- 
eters of the cell (number of MTs) and the virus (binding 
and unbinding rates, the mean velocity V of the directed 
motion and the diffusion constant D). 
General Methodology. To derive an expression for b(x), 
we consider the motion of a virus between the moment it 
enters the cell at the outer membrane and the moment it 
reaches the absorbing boundary dN a . Its motion alter- 
nates between free diffusion, for a random time t, until 
it hits a MT and binds. It continues in a directed motion 
along the MT for a mean time t m and a mean distance 
d m = \\V\\t m , until it is released and resumes free diffu- 
sion. The steady state drift is chosen to be constant for a 
sufficiently small step, such that the mean time T + t rn to 
the first release at a point Xf is the same as that predicted 
by @ (see FIG. [I}. This approach leads to explicit ex- 
pressions for the steady state drift for two-dimensional 
radial and cylindrical geometries. 

The steady state drift for a two-dimensional radial cell. 
We consider a two-dimensional cell cytoplasm which is 
an annulus Q of outer radius R and inner radius 5 (nu- 
clear surface) with N MTs radially uniformly distributed. 
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FIG. 2: Two dimensional radial cell with radially equidis- 
tributed MTs. We show a virus trajectory alternating be- 
tween bound and diffusive periods in cytoplasm. 
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FIG. 3: A fundamental step in Q. The virus starts at a 
radius ro, with an angle uniformly distributed in [0;G], it 
diffuses freely during a time r(ro) until it binds to a MT at 
a mean radius r(ro); it has then a directed motion over a 
distance d m = ||V||t m before being released randomly at a 
final radius r/. Mean times of each piece of the fundamental 
step are written inside parenthesis. 



They irradiate from the nucleus to the external mem- 
brane and the angle between two neighboring ones is 
= The two-dimensional approximation applies for 
culture cells which are flat [3] due to the adhesion to the 
substrate. In that case the thickness can be neglected in 
the computation. Before reaching a small nuclear pore, a 
virus has an intermittent dynamics, alternating between 
diffusing and bound periods (see FIG.[2]). Because the 
MTs are uniformly distributed, we consider the funda- 
mental domain Q defined as the two-dimensional slice of 
angle between two neighboring ones. In f2, the fun- 
damental step described above is as follows: the virus 
starts at a radius ro with an angle uniformly distributed 
in [0; 0], it binds to a MT at a time r(ro) and at a radius 
f(ro). On the MT, it has a radially directed movement 
towards the nucleus during a mean time t m and over 
a distance d m = ||V||t TO . Finally, the virus is released 
with a 0-uniformly distributed angle at a final radius 
77 = f(rrj) — ||V||i m (sec FIG. [3]). In most eukaryotic 
cell large asters, there are from 600 to 1000 MTs Q. 
We can estimate the average number N of MTs per cell 
cross section as follow: for a cell thickness h w 9/zm, 
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Q, an interaction range 7 « 50nm between the MTs 
and the molecular motors 0], and for the AAV diameter 
d = 30nm Q, we obtain for a radial MT organization 
in a thin cylindrical cell, that the range of N is between 
600(2 J +d) and 1000( ^ +rf) , that is 9 to 15. We are thus 
in a regime where << 1. For ro < R, by neglecting 
the reflecting external boundary at r = R, becomes an 
open wedge and thus using the standard methods from 
15, 16], we obtain 



@2 Q2 
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In radial geometry, b(x) = 6(r)y^yj and the MFPT u(ro) 
of a virus starting at ro and ending at position r/, de- 
scribed by equation Q satisfies [l7j : 

(7) 
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where we approximated 6(r 
equation {2} is 



u(^o) 



L»Aw - 6(r )Vu 
(P) = and 1/(77) = °> 

by b(ro). The solution of 
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For D << 1, using the Laplace method, 
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(9) 



To obtain 



Thus, in first approximation, u{r ) w "gzpff ■ 
the value 6(ro), we equal the MFPT u(ro) from ro to re- 
computed from equation |(2J) with the one obtained from 
an intermittent dynamic: r(ro) + t m . Consequently, we 
get: 



b(r ) 



r_p - Tj _ a 
r(r ) + t m t m + r, 
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(10) 



Tests against Brownian simulations. We impose re- 
flecting boundaries at the external membrane r = R 
and we tested the theoretical steady state distribution 
against the one obtained by running empirical intermit- 
tent Brownian trajectories in the pie wedge domain. For 
a potential field, the steady state distribution satisfies 
DAp — V[bp] = in f2 with reflecting boundary condi- 
tion J(x, t).n x = on dfl. The distribution p in a two- 
dimensional radial geometry is: 



p{r) 



e d 2nrdr 



(11) 



which should be compared to the distribution of [zj . The 
potential $ of b = — is obtained by integrating equa- 
tion (flTjj) with respect to r, 



S(r) 
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(12) 
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FIG. 4: Steady state distributions. Dashed line: virus distri- 
bution with the effective drift b{r) (|10[) ; solid line: Empir- 
ical steady state distribution obtained by running 10, 000 in- 
termittent Brownian trajectories. The cell radius is R = 20/im 
and9 = £. 



In FIG. [H we plotted the steady state distribution given 
in (jlip against the distribution obtained by the intermit- 
tent empirical equation {TJ. The parameters are chosen 
such that the viruses move towards the nucleus (observed 
in vitro, loaded dynein moves during Is over a distance 
of 0.7/xm [3), we thus take t m — Is and d m — 0.7/um; 
furthermore, the diffusion constant is D = 1.3/xm 2 s _1 as 
observed for the Associated- Adeno- Virus Q. The nice 
agreement of both curves, which is the central result 
of this communication, confirms that our coarse grained 
method accounts well for the switch system (UJ. 
Computation of Pn and tn ■ We derive now asymptotic 
expressions in the small diffusion limit D << 1, for the 
probability P/v and the mean time rjy a virus arrives at a 
small nuclear pore. We apply Laplace's method in formu- 
las for a radial geometry. When the degradation rate 
k{r) is taken constant, equal to fco in the neighborhood 
of the nucleus r = S and when 12d m > rO 2 , b{r) > so 
that $ reaches its minimum at r = 6, we get 



Pn = 



b(S) 



and tn 



ln{\) 28 



ln(\)25k Q + b{5) In (\) 26k + b(S)' 

A Taylor expansion for Q <C 1 gives that 



Pn 



tn 



K 



1 - 



k {dm + K) 



KS {d m S + Dt m ) 
l2Dt m d m {d m + K) 
5{d m 5 + Dt m ) 



\2Dtm {d m + K 



6 2 ) (13) 
9 2 \ (14) 



where K = 2koSt m In (-). 

We can now propose the following predictions: because 
nuclear pores occupy a fraction e = 2% [l9[ of the nucleus 
surface (radius S = 8«m) and the measured degradation 
rate for plasmids [20( is k = g^s -1 , we obtain from 
formula [TSril that 



P w ps 94.3%, t n « 205s. 



(15) 



We conclude that the infection efficiency is very high, 
while the mean time to reach a nuclear pore is of the 
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FIG. 5: Dendrite cross-section. The N MTs are thin cylinders 
uniformly distributed inside the dendrite. 



order of 3 minutes. It is interesting to compare this time 
with the 15 minutes reported in Q, which accounts for all 
the viral infection steps from the entry to the final nuclear 
import. This difference between the two times indicates 
that the phase where the virus is inside an early endosome 
(EE) may last 10 minutes. Indeed, the endosomal phase 
ends once the EE has matured into a late endosome (LE) 
[2lj . which lasts approximately 10 minutes 22}. To finish, 
we shall note that a free diffusin g; v irus would reach a 
nuclear pore in about 15 minutes [111 ]. 
The cylindrical geometry. Many transports mechanisms 
such as viral (herpes virus |23|) and vesicular occur in 
long axons or dendrites, which can be approximated as 
thin cylinders (radius R and length L). To derive a quan- 
titative analysis of viral infection in that case, we fol- 
low the method described above and compute the steady 
state drift that accounts for the directed motion along 
MTs. We model the N MTs parallel to the dendrite prin- 
cipal axis as cylinders (radius e << R, Length L). The 
cross-section 51 of the dendrite is shown in FIG.O Due to 
the cylindrical symmetry, for any position x, the steady 
state drift b(x) is equal to Bz where B is a constant and 
z the principal axis unit vector along the dendrite. In 
a small diffusion approximation, the leading order term 
of B is equal to the effective velocity [!, H| : B = t d ™ T , 
where t m is the mean time the virus binds to a MT, 



dm = 1 1 V| \t m the mean distance of a run and r the MFPT 
to a MT. To derive an expression for r, we consider the 
cross-section f2 and impose reflecting boundary condi- 
tion at the external membrane of the dendrite (r = R) 
and absorbing ones at the MTs surfaces. In long time 
approximation, for a MTs radius e << 1, r is asymptot- 
ically equal to where A is the first eigenvalue of the 
Laplace operator in f2 with the boundaries conditions de- 
scribed above (fl7| p. 175). The leading order term of A 
as a function of e is [IJ] A = j^j^rn » where |fi| = kR 2 . 

Thus, the MFPT to a MT is r = ^ = ^nd < and the 
steady state drift amplitude B is given by 



B = 



2NDd n 



2NDt r . 



R 2 ln 



(16) 



We conclude that in the limit t m <C r, the effective veloc- 
ity is proportional to the number of MTs: B w N \ , 

as already observed in Q. 

Conclusion. Intermittent dynamics with alternative peri- 
ods of free diffusion and directed motion along MTs char- 
acterizes many cellular transports. We have developed a 
model to estimate a steady state drift such that the in- 
termittent dynamic can be described by an over-damped 
limit of the Langevin equation. Our method gives explicit 
results in two-dimensional radial cell and in a cylindri- 
cal dendrite or axon. The steady state description of the 
movement enables us to estimate the probability a virus 
reaches alive a small nuclear pore and its mean time. Be- 
cause viruses are very efficient DNA carriers, understand- 
ing and quantifying their movement in the cell cytoplasm 
would be very helpful for designing synthetic vectors [25[ . 
In a future work, it would be interesting to derive steady 
state drifts for three dimensional geometries. 
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